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ABSTRACT 



Context. Sun-grazing comets almost never re-emerge, but their sublimative destruction near the sun has only recently been observed 
directly, while chromospheric impacts have not yet been seen, nor impact theory developed. 

Aims. We seek simple analytic models of comet destruction processes near the sun, to enable estimation of observable signature 
dependence on original incident mass M„ and perihelion distance q. 

Methods. Simple analytic solutions are found for M(r) versus q and distance r for insolation sublimation and, for the first time, for 
impact ablation and explosion. 

Results. Sun-grazers are found to fall into three (Mg, q) regimes: sublimation-, ablation-, and explosion-dominated. Most sun-grazers 
have Mg too small (< 10" g) or q too large (> \.Q\Rq) to reach atmospheric densities (« > 2.5 x 10"/cm^) where ablation ex- 
ceeds sublimation. Our analytic results for sublimation are similar to numerical models. For q < \.Q\Rq„M„ > 10" g, ablation 
initially dominates but results are sensitive to nucleus strength P^ = 10* Pg dyne/cm^ and entry angle to the vertical. Nuclei with 
Mo < lO'^CPe sec (jyf g are fully ablated before exploding, though the hot wake itself explodes. For most sun-impactors sec (^ » 1 
(since q ~ r,), so for q very close to r, the ablation regime applies to moderate M,, ~ lO'^^'^fg g impactors unless P(, < 0.1. For 
higher masses, or smaller q, nuclei reach densities n > 2.5 x IQ^^PJcm^ where ram pressure causes catastrophic explosion. 
Conclusions. Analytic descriptions define (Mo, q) regimes where sublimation, ablation and explosion dominate sun-grazer/-impactor 
destruction. For q < l.Ol^g, M^ > 10" g nuclei are destroyed by ablation or explosion (depending on M„ cos' ip/Pc) in the chro- 
mosphere, producing flare-like events with cometary abundance spectra. For all plausible Mo,q and physical parameters, nuclei are 
destroyed above the photosphere. 
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1. Introduction 

Close sun-grazing comets are discovered almost daily by 
white light coronographs (e.g. SoHO LASCO). Most are 
small and fully sublimated by insolation at a few solar 
radii Rq while almost none have re-emerged (Marsden 
2005). The majority have perihelion distances q well above 
Rq but some have q ^ or < Rq (e.g. Marsden 2005 and 



http://www.ininorplanetcenter.org/mpec/RecentMPECs.htnill) 



The death of a comet at r ~ ^0 has been seen directly only very 
recently (Schrijver et al 201 1) using the SDO AIA XUV instru- 
ment. This recorded sublimative destruction of Comet C/2011 
N3 as it crossed the solar disk very near perihelion ^=1.1 39Rq. 
The next challenge in studying the demise of close sun-grazers 
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will be to catch one of the even rarer cases of chromospheric 
impact iq < \.Q\Rq, Mo > lO" g). As anticipated by Weissman 
(1983) and shown below, these undergo explosive destruction 
in the dense chromosphere. Understanding the destruction 
processes, and their radiation signatures, are essential steps in 
searching for and modeling these. 

The processes leading to sublimation of the icy conglomer- 
ate mix (Whipple 1950) of cometary nuclei, and in some cases 
their splitting and fragmentation, were considered by Huebner 
(1967), Weissman and Kieffer (1981), Weissman (1983), Iseli 
et al. (2002), Sekanina (2003) and others. These models essen- 
tially solve for the insolative sublimation mass loss rates of icy- 
conglomerate mixes (the dust being caiTied away in the flow of 
these evolved components). They variously allow for the com- 
plicating factors of rotation, albedo, insulating surface dust lay- 
ers, radiative cooling, interior thermal conduction, and fragmen- 
tation by tidal, thermal and volatile explosion effects. Huebner 
(1967) and Iseli et al. (2002), for example, found that, for high 
sublimation rates near the sun, these effects were secondary, 
the mass loss being reasonably approximated by a pure subli- 
mation description : mass loss rate - heating rate/latent heat. 
Using SoHO data, Sekanina (2003) addressed in detail how mass 
loss rates and fragmentation relate to cometary light curves via 
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atomic line emission (e.g. by sodium) and by dust scattering of 
sunlight, though emphasising that most of the mass remains in a 
primary fragment. None of these studies considered ablation or 
ram-pressure driven explosion due to solar atmospheric impact 
(which we show below are negligible till r < \.Q\Rq) though 
Weissman (1983) had remarked "ultimate destruction of the nu- 
cleus [of sun-impacting comets] likely results from the shock of 
hitting the denser regions of the solar atmosphere" . It is well 
known in the planetary physics community that ablation and ex- 
plosion are central processes in comet-planetary atmosphere im- 
pacts (e.g. Carslon et al. 1997) 

Here we revisit the theory of sun-grazer sublimation then 
develop the first estimates of the much higher rates of mass 
loss by ablation for the rarer cases which reach < \.Q\Rq. 
This regime (Sections 3.4, 6.1, 6.2) resembles that of comet- 
planetary impacts though with some differences. One is that al- 
most all sun-grazing comets belong to secondary comet groups 
(mostly Kreutz) formed from primordial comets. A significant 
number of the latter must have q < Rq (Hughes 2001) and large 
enough mass (Mo > 10" g - see Section 5, Eqn. (IZTT i) to sur- 
vive sublimation down to the intense ablation/explosion regime. 
However only a small fraction of group comets (Biesecker et 
al. 2003, Knight et al 2010) come that close or are that mas- 
sive, larger sun-grazers mostly having q > 1.5Rq. Comet C/2010 
E6 (STEREO) discussed by Raftery et al. (2010) came close, 
having q a: 1.02Rq while the destruction of Comet C/2011 N3 
seen by SDO (Schrijver et al 201 1) on the solar disk was solely 
by sublimation as it had q * I.IARq. Most of the few group 
comets reaching the ablation/explosion regimes will have small 
mass and very shallow incidence angles (since q ^ Rq) with im- 
plications for the relative importance of ablation and explosion. 
By analogy with models of planetary atmospheric impacts (e.g. 
Carlson et al. 1997) , chromospheric impactors are expected to 
lose mass, momentum and kinetic energy very rapidly, the soar- 
ing temperature and pressure leading to destructive detonation 
of the nucleus and its wake. This 'air burst' should be followed 
by explosive expansion of the resulting 'fireball' of cometary 
and solar atmospheric material. The total masses and energies 
of these explosive chromospheric events are in roughly the same 
range as those of solar magnetic flares so they could be termed 
'cometary flares'. 

Section 2 discusses relevant nucleus parameters values while 
Section 3 estimates the relative importance and time/distance- 
scales of nucleus mass loss, explosion, and deceleration pro- 
cesses. We then consider in some detail approximate ana- 
lytic treatment of mass loss (Section 4) purely by sublimation 
(Section 5) and by ablation and explosion (Section 6). 

2. Comet Nucleus Parameters 

2. 1 . Masses and Densities 

In the past there was wide range of opinions over the masses 
Mo and mean mass densities p of the icy-conglomerate mix 
thought to comprise comet nuclei, and on how these quanti- 
ties varied between different comets (Whipple 1950). We sim- 
ply take the incident mass M,, to be a parameter ranging from 
the faintest comets detected almost daily to the largest super- 
massive comets. Concerning the density p, MacLow and Zahnle 
(1994) said the density of comets is the stuff of speculation [prob- 
ably] 0.3 - 1.0 glcrr? but with extreme claims of 0.01-5. In his 
sublimation modeling Sekanina (2003) used p - 0.15 gjcvD? . 
More recently, however, Weissman and Lowry (2008) critically 
discussed results from a range of methods and data, including 



Deep Impact. All results have a considerable error range, but 
they conclude that p - 0.6 ± 0.2 g/cm^ is the most likely range 
for the comets they studied (cf Richardson et al. 2007). Here we 
therefore also take p as a parameter, and express it relative to 
that of (water) ice pice - 0.9 glcvc? via the dimensionless scaling 
parameter 



P = PiPic 



(1) 



In discussing results we mainly use a mean value p - 0.5 (p ^ 
0.45 gjcnr') though showing how results, such as for Mo, can be 
scaled for other values (see Eqn. ( l20b and Figure 1). 

2.2. Latent heat of sublimation 

The major nucleus parameter determining its rate of mass loss 
for power input P (erg/s) is the mean latent heat X (erg/g) of sub- 
limation/ablation for the cometary ice-conglomerate mix. The 
value for water ice is about 2.6 x 10"^ erg/g. This is of order the 
relevant value when the mass lost is in the form of water vapor 
whose pressure sweeps dust and other non-volatiles away with 
it. However, in situations of intense ablation and explosion, all 
components are eventually vaporized and the relevant value is 
the weighted mean over all mass components, including volatile 
and refractory ones. Sekanina (2003) finds values for what he 
terms 'effective latent heat of erosion' tobearound0.3of the wa- 
ter ice value. However Huebner et al's. (2006) mass composition 
[Silicates : organics : carbonaceous: ices ] - [0.26:0.23:0.086: 
0.426] suggests refractory silicates to be the most important non- 
ice component. Mendis (1975) give X - 2.30 x 10'" erg/g for 
silicates, which is close to that of water ice, while organics will 
lower the mean slightly. Chyba et al (1993) adopted 2.3 x 10'° 
erg/g for comets, 5 x 10'" erg/g for carbonaceous and 8 x 10'" 
erg/g for stone/iron bodies and just 10" erg/g even for solid iron. 
Here we use the water ice value for X as a reasonable first ap- 
proximation, but include a dimensionless scaling parameter 



-C = Z/£ic 



(2) 



to allow application to other X values though adopting X ~ 
1, (X ~ X/cf), in most numerical evaluations. To convert to other 
values see Figure 1 and Eqn. (l20t . 

2.3. Physicai Strength and Sound Speed 

The effects on the nucleus of the pressure associated with the 
processes driving sublimation/ablation depend on how strong 
the nucleus material is and how fast pressure waves propagate 
through its volume. In the initial stages the relevant strength is 
that of the loose low density ice-conglomerate mix, estimated 
by Zahnle and MacLow (1994) to be a pressure (^ energy den- 
sity) of Po ~ 10* dyne/cm^(= erg/cm^) while values discussed 
by Chyba et al (1993) range from 10^ dyne/cm^ for very solid 
asteroid-like bodies to 10^ dyne/cm^ for fluffy snowy structures. 
The last value is also that estimated by Richardson et al (2007) 
for Comet 9P/Tempel 1 using Deep Impact data, while tidal 
disruption data suggest the possibility of even lower values for 
some comets. Here we allow consideration of a range of values 
by parametrically setting P - 1Q^P(, dyne/cm^ and specifically 
discuss the range Pf, = 0.1 - 1. For the sound speed, Zahnle and 
MacLow (1994) suggest c., w 1 - 3 x 10^ cm/s for the initial 
state nucleus material. Values for a range of more or less solid 
materials do not differ greatly. 
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2.4. Nucleus shape and 'size' 

For a given density the mass loss, deceleration and explosion 
rates for a nucleus depend on its volume V and effective area 
A which depend on its size scale 'a' and shape. Since nuclei 
have irregular and uncertain shapes, as evidenced by images and 
light curves, rather than treat precisely the hypothetical case of a 
sphere, we will use A ^ a^,V ^ a^ for an effective mean size a. 
This formulation gives A and V correctly to within 10% for the 
case of a sphere of radius b if we adopt a ^ 1 .67b. We use this 
later when referring to spherical cases for illustration (Sections 
5.2.2 and 5.3). 



3. Rough Comparison of Nucleus Mass Loss, 
Deceleration and Explosion Processes 

3.1. 'Evaporative' Mass Loss versus Deceieration 

Whether mass loss from the nucleus is by sublimation or by ab- 
lation, the total energy needed to vaporize the whole nucleus is 
£>van - Mo£. while the energy needed to stop it is £i,„ ^ MoV^/l. 



-'vap 



The ratio is tiny, &vapl&km ~ '^■CIvq = 4 x 10 ■'-C for the solar 
escape speed 



V = vo = (IGMJRs) 



1/2 



618km/s 



(3) 



with G,Mq,Rq the gravitational constant, solar mass, and solar 
radius respectively. Thus total vaporization can occur well be- 
fore the nucleus decelerates significantly, though ram pressure 
can result in explosive destruction dominating. Equivalently, to 
conserve momentum in slowing down, the nucleus must en- 
counter an atmospheric mass comparable to M„. In doing so it 
absorbs far more energy than needed to vaporize it. The mass 
per unit area (g cm~^) of a comet nucleus is 2^ ^ Mg/a^ - 
mI'^P^I^ = lQ^{Moll0^y^{pfl^ g cm-2 while that of the sun's 
atmosphere even down to the photosphere is only 2© a; 1 g cm"^. 
Consequently, unless they explode, increasing the deceleration, 
only objects of < lO"* g or so would be much decelerated by the 
mass of the sun's atmosphere down to the photosphere though, in 
practice, they would be vaporized earlier by sublimation and/or 
ablation. Note also that, at infall speed Vq, in the frame of the 
nucleus a solar atmospheric proton has kinetic energy ~ 2 keV 
which is enough to knock off about 4000 water molecules, or 
about 10^ times its own mass. 

This situation of initial ablative mass loss without significant 
deceleration is paralleled by the dynamics in planetary atmo- 
spheres of meteors - e.g. McKinley (1961), Kaiser (1962) - and 
at least of the initial (high altitude) stages of comet-planet im- 
pacts - e.g. Shoemaker Levy 9 with Jupiter - cf. Section 6. 

In the sublimation regime the radiation has negligible effect 
on the nucleus speed, even though it delivers a large amount 
of energy. Specifically the ratio of the radiation pressure force 
Frad to the gravitational is (ignoring the correction factor 1.0- 
2.0 for albedo) y^ < ^'"^-- . The radiation pressure F,ad/c - 2 

dyne/cm^ is also tiny compared to the nucleus strength P^ so 
causes no direct explosion effect. It is also straightforward to 
show that, during perihelion passage, the rocket effect of mass 
loss leaving the nucleus anisotropically has negligible effect on 
V during vaporization. This is because the very high nucleus p 
value implies very low mass loss speed u and hence momentum 
flux. Consequently, the velocity of the nucleus during its vapor- 
ization is well approximated by that of a gravitational parabolic 



orbit, viz., in orbital plane polar coordinate (r, ff) and for perihe- 
lion distance q 



(Rr.\"', 



Vf — r - 






1/2 



(4) 



3.2. Deceieration of subiimated/abiated mass 



While, throughout the total vaporization lifetime of the nucleus, 
its velocity is well described by Eqns. (|4]i, this is not true of 
the material it loses. Small particles leaving the nucleus do not 
obey Eqns. dUi, non-gravitational accelerations on them being 
very important. These include F,-ad on sub-micron particles and 
atmospheric drag on atoms and ions. The Coulomb collisional 
stopping column density for a proton of speed Vq « 618 km 
sec"' (~ 2 keV) in a hydrogen plasma is around Ns - lO'*" cm"^ 
(Emslie 1978). Thus the stopping distance at number density n 
is li(cm) a N/n » lO'^'/n, i.e. just 1 km in the corona and 0.1 
mm in the photosphere. So ablated dust and ions stop abruptly 
and form an exploding wake as they blend with and heat the 
atmosphere, creating large local enhancements of heavy element 
abundances. 



3.3. insolation Sublimation case ^ > r,, Mo < 10" g 

Very near the sun, the timescale for sublimation of the whole 
mass is T^ut > MoZ/alTo ~ 4 x lO^A/Jf p^/a^ ^^ ^j^^j.^ ^^ ^ 
6x10'" erg/cm^s is the bolometric photospheric energy flux and 
M\2 = Mo/lO'^ g. The corresponding distance scale is roughly 
dsiib ~ VqTsui, ~ 3.5Rq X Mj2 p^^^X. cm as first seen directly 
by Schrijver at al. (2011). Total sublimation occurs in a close 
perihelion passage roughly for dsui, ^ ^0 which is the case for 
masses M^ < 3 x lO^'/p^X^^ g or about 10"g forp = 0.5, X = 1. 
This shows the majority of close sun-grazers {q > I.QIRq) to 
have Mo < 10" g as they never re-emerge. 



3.4. Impact ablation and ram pressure explosion cases 

3.4.1. Background 

The problems of comet, asteroid, and even large meteoroid im- 
pact with planetary atmospheres closely parallel those of a so- 
lar impact, though the parameter regimes are rather different, 
and there is no planetary equivalent of insolation. Unlike the so- 
lar case, planetary impacts have been addressed in detail, early 
work including that of Revelle (1979) and others. Progress was 
greatly accelerated in anticipation, and in the aftermath, of the 
collisions of fragments of Comet Shoemaker-Levy 9 with Jupiter 
in July 1994. Many authors (e.g. Chyba et al 1993, ChevaHer 
and Sarazin 1994, Zahnle and MacLow 1994) developed semi- 
analytic and numerical models to predict what should be ex- 
pected of these impacts. Others, notably MacLow and Zahnle 
(1994), Field and Ferrara (1995) and Carlson et al. (1995, 1997), 
developed models further by drawing on actual event data. In 
particular the Carlson et al. (1995, 1997) 'heuristic model' com- 
bined data analysis, numerical simulations, and observational 
inputs. They identified and addressed both the 'bolide initial 
phase' (which we have called the detonation or airburst) and 
the subsequent 'fireball' exploding out and up from the nucleus 
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and its wake. Detailed observations of the G-fragment impact 
were synthesised with their numerical simulations of the fireball 
(e.g. their Figure 7). Chyba et al. (1993), Zahnle and MacLow 
(1994) and MacLow and Zahnle (1994) showed that the airburst 
phase is dominated by ablation only in the upper atmosphere, 
and below that by catastrophic pressure-driven disruption. In the 
Jupiter case, the latter is true for all impacting masses except 
those very much smaller than that (^ 10'^ g) of the Shoemaker- 
Levy 9 fragments in 1994. We summarize here the essence of 
the arguments leading to these conclusions so that we can adapt 
them to the somewhat different conditions of solar impacts. 

3.4.2. Rough comparison of the importance of impact 
processes 

The relative importance of ablation, explosion and deceleration 
in the initial destruction of the nucleus depends on atmospheric 
scale height H as well as on the nucleus mass, speed and entry 
angle. They can be expressed in terms of the atmospheric col- 
umn densities NicmT^) - J nds traversed along the path s of 
the nucleus after which each processes is complete. (The vertical 
column density N - N cos (p for constant path angle (p to the ver- 
tical) . First, to deliver enough energy to ablate the whole mass 
(without expansion) requires an encounter with column density 
Nahi such that Nabi/J^nipV^/l ^ HMo/al where irip is the proton 
mass and fi the mean mass present per H (atoms and ions) in 
units of rrip. This implies 

Khiicm-^) = Nm sec0 ^ £M'fp^'^/(jjmpvl/2) 

a; 6 X lO^^Mj'f Xp2/3 (5) 

Total deceleration of the nucleus (without expansion) requires 
Ndec such that Ndecf^mp ~ Mo/a^ implying 



Ndecicm ^) 



I^Mp 



: 4 X lO^^Mjf p2^^ 



,27 ,.1/3 



3 X WM[i; (6) 



which is » Nabi- For internal strength (energy density) Pc - 
lO^f 6 dyne/cm^, nuclei will undergo hydrodynamic flow once 
they reach the depth where atmospheric ram pressure fiiiipV^/l 
approaches the value Pc- This implies an atmospheric proton 
density n > tipres = «** ~ ^Pc/iumpV^j. The corresponding col- 
umn density Npres along the path for atmospheric scale height 
H = lO^Hs cm is 



,-2' 



NpreAcm ) = HpresHseC 



2PcH sec (I 

2AxlO^'^P6Hssec<p 



(7) 



The ratio of critical N for ablation end to ram pressure onset is 
(Eqns. ailli) 



^„, 



0.4P6^8 sec 4> 0.7 PeHi sec c, 



K 



abl 



M''£p^l^ 



M, 



1/3 
12 



(8) 



where the last expression is forp = 0.5, X - 1. There is thus a 
maximum mass M™j'((p) above which destruction of the nucleus 
is driven by explosion rather than ablation, namely (setting ratio 
(O to unity) 



M'Z'id),) = M„ sec^^ * 6 X 10'° 



^abl 



P6Hgsec4 



Xp2/3 

2.4xlO'i[P6//8sec0]3g 



the last value again for p - 0.5, X. - I- (Note - In the rest of 
this paper we use subscripts * and „ respectively for values of 
quantities at the points where ablation first exceeds sublimation, 
and »» where explosion first becomes dominant.) 

Thus although ram pressure sets in at a lower n for the sun 
than for Jupiter (because v is higher) the same v dependence oc- 
curs in the ablation expression so the ratio ([8]) just scales as 1 /H 
and reaches unity only at a higher n in the solar case. For solar 



//g ~ 0.5 we find n 



pres 



end 
ahl 



n„ = 2.5 x W^Pk cm-^' while 



(9) 



M„ * 3 X 10'°P^ g and M'2' ~ 3 x 10'°(P6 sec 4>f g. In the case 
of impacts with Jupiter //g ~ 0.01 - 0.1 so, for relatively steep 
impacts (sec 4> < 2), Eqns. ([8] |9]l show that ram pressure will 
cause explosion before ablation is complete for Mo > 2 x lO^f ^ 
g. Thus, as argued by Chyba et al (1993), Zahnle and MacLow 
(1994) and Carlson et al. (1997), for steep planetary impacts the 
ablation process is secondary to pressure-driven explosion ex- 
cept for small masses (or extremely shallow entry) and during 
the initial entry phase. In the case of solar chromospheric im- 
pacts, however, i/g - 0.5 is larger so that the maximum ablation- 
dominated mass (after sublimation losses) is around 3 x 10'"Pg 
g for steep entry. However, most sun-impactors have q k R^ and 
sec(^ » 10. Thus the maximum ablation-dominated mass for 
these is » 10'"*^^ g as discussed further in Section 6.1. 

In the planetary impact literature two factors are discussed 
which act to limit ablation and we show briefly here why they are 
less relevant to the solar case. The first is possible overestima- 
tion of the drag coefficient - we have taken the nucleus effective 
area to equal its geometrical area, which is an overestimation for 
fluid flow round an obstacle. However, for the solar case the nu- 
cleus speed is so high and the atmospheric density so low that 
we can reasonably use a kinetic description and treat the parti- 
cles as impacting the nucleus directly rather than flowing around 
it. We note, however, that our kinetic approach may break down 
toward the lower chromosphere relevant to very heavy nuclei 
and steeper incidence angles. There the high density results in 
a shorter collisional mean free path and more fluid like behav- 
ior with only indirect ablation via the radiation from a stand- 
off shock. This can decrease the drag coefficient considerably so 
that our ablation and deceleration rates are overestimates there. 
Secondly, it has been argued in the planetary case (e.g. Zahnle 
and MacLow 1994) that ablation is limited by the loss of en- 
ergy needed to ionize hydrogen. However, the solar atmosphere 
is around 100 times hotter than Jupiter's atmosphere and is al- 
ready substantially ionized at least in the upper chromopshere. 
Further, solar impacting protons have keV kinetic energies - far 
higher than H ionization energy (13.6 eV) and than the 10 eV 
value for Jupiter so direct impact ionization is very effective in 
the sun but not Jupiter. Again, however, heavier and steeper entry 
nuclei reaching the deep chromosphere are in a cooler medium 
of lower ionization so energy lost to ionization may become sig- 
nificant and our ablation rate an overestimate. 



4. Theory of Nucleus Sublimation/Ablation 

4.1. General Formulation 

Having discussed orders of magnitude we now want to model 
the evolution of evaporative mass loss along the comet path. We 
take the nucleus to have constant p and to shrink due to mass 
loss, neglecting for now any pressure driven explosion or dis- 
tortion which we address later. The changing M{r) and a{r) are 
then functions of distance r from the solar centre. We treat mass 
loss (by sublimation or ablation) via the simple approximation 
that a heating/energy deposition rate of y' erg s ' results in mass 
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loss rate M g s"' = PIJl with the relevant p,£. values aver- 
ages over the icy conglomerate components (Sections 2.2.1 and 
2.2.3). For the sublimation case the dominant power required for 
loss of mass from the nucleus is that for sublimation of the ices 
only, the solid dust particles released escaping with the vapor 
outflow. That is, fairly near the sun at least, the power going into 
sublimation is larger than that going into further heating off'set 
by cooling - e.g. Weissman (1993), Sekanina (2003). (Far from 
the sun our neglect of cooling etc will overestimate the mass 
loss rate). For the ablation case, heating is intense and vapor- 
izes all components, including the most non-volatile ones but 
(Section 2.3), the overall X is similar to X/cf We further assume 
that all mass lost clears the nucleus quickly enough for incident 
energy fluxes to act directly on it - i.e. the mass loss does not 
obstruct the mass loss driver. To justify this, we note that la- 
tent heat £. is equivalent to a speed ~ 2X'^^ ~ 2.3 x 10^ cm/s. 
If vaporized material outflow occurs at this speed then, on the 
timescale ///v© cos ~ 10 s relevant to ablation, it clears nuclei 
of a,, < 2.3 X 10^ or M,, < 10'^ g so our assumption is justified 
except for very large masses. These are in any case destroyed 
by explosion rather than ablation. In the case of sublimation the 
mass loss rate is much lower (Section 3.3) so the approximation 
is very good. 

Then, for an incident energy flux T (erg cm^^s '), V ~ a^T 
and, with a^ - (M/p)^^^, and constant p and shape, we get 



M ■ 



dM 
dt 









(10) 



Note that we consider only mass loss in the one hemisphere 
(cross section w a^) containing flux T since rotation times are 
long (Samarachinha et al. 2004) compared to the short near-sun 
mass-loss time scale. Using Eqn. (JUi to express dt - -drjvr we 
can eliminate t from Eqn. ( fTOb and express it as a differential 
equation for the fractional mass m - M{r)IMo remaining as a 
function of r, viz. 



1 dm R] 



3/2 



1 



'F(r)r 



1/2 



^2/3 dr vo Ml'^£pyHl-q/ry'^ 



(11) 



This can be integrated for any specified form of Tir) with 
boundary condition m(r — > oo) = 1. We have formulated the 
problem in terms of r dependence since this allows analytic treat- 
ment. Conversion to the time domain can be done numerically 
via Barker's equation (cf. results in Porter 2007, 2008). 

4.2. Heating Energy Flux Terms and regimes 
- (a) absorbed insolation energy flux for albedo «: 1 is 



Tr 



rad 



fir) 



Anr^ Anr^ 



5xlO»'(^) ergcm-2s-i(12) 



where f{r) is the factor by which the heating rate at r exceeds 
our estimate in Eqn. (fT2l) where we use the approximation of 
near radial insolation. As r — > /?© non-radial insolation from 
the finite solar disk increases and for a sphere just above the 
photosphere the maximum is /(/?©) = 2 (= 4 below the pho- 
tosphere) but this can increase the mass required to survive 
insolation 8-fold, fir) declines above r - Rq, is shape de- 
pendent, and partly offset by the cooling processes we have 
neglected but our use of / ~ 1 must be borne in mind in 
considering our results. 



(b) ablation energy flux Impacting solar protons have a col- 
lisional mean free path (mfp) in nucleus material (n ^ 10^^ 
cm"-') of order lO"*" cm (cf. Brown 1972, Emslie 1978) while 
in the atmosphere at n > n, = 2.5 x 10", the mfp is 
~ 4 X 10^ cm, much bigger than most nuclei. Furthermore, 
the nucleus speed is highly supersonic (~ Mach 50 in the 
chromosphere). Thus the bombarding atmospheric protons 
have neither enough time nor undergo enough collisions to 
behave as a fluid. We therefore treat the problem as one of 
surface ablation by this particle kinetic flux rather than as a 
fluid flow. The particle bombardment energy flux is then 



1 .. 3, ^ 1 ,, 3(Ref^ 



Tcoiiir) = -P&(y)v (/) = -/UOTp 



= 2.7x10 



riQ \ r I 



3/2 



erg cm ^ s ' 



(13) 



(cf. meteor mass loss in the earth's atmosphere - e.g. Kaiser 
1961, Chapter 8). Here n{r) is the total (atom + ion) number 
density (cm^^) of hydrogen, nip the proton mass and fi ^ 
1.3 the mean total mass present (including helium etc) per 
hydrogen (so thatp© = fJ.nmp) while n© x 10'^ cm"-' is the 
value of n at the photosphere (r = Rq) which we take as 
reference point. 
ratio of ablation to sublimation energy fluxes is 

Z£^.4xlO^(4''^.(4^'^ (14) 



n« 



where n* = n(r,) ^ 2.5 x lO^cm""* is the n value at r = 
r» = I.QIRq where Trad,'Fcoii are equal. At the photosphere 
in - Hq- 10'^ cm"-'), Tcoii ~ ^O^f^rad which is equivalent to 
a black body energy flux at T > 100, 000 K ! 

From Eqn. ( fT4l l we see that sublimation dominates over ab- 
lation only at distances where n{r) < n^ ^ 2.5 x 10" cm"-' 
where they are equal. This value n - n, is in the (upper) chro- 
mosphere where the density rises very steeply with depth with a 
fairly constant temperature (of order 10,000 K) and hence mean 
exponential density scale height H ^ 5 x IQ^ cm w 7 x lO""*/?©, 
viz. 



n(r) 



„-(r-Ro)/H 



(15) 



In this approximation, the values r - r.,,n - n^ occur at a height 
z above the photosphere of z - Zt - r^ - Rq ^ //ln(4 x 10^) ~ 
6500km « 0.01/?o. Since /i = HIRq ^ 7xl0"^overavery smaU 
r range (a few H) the dominant energy flux changes rapidly from 
radiative at z - z* » // to collisions at z - z, <« -H, a z range 
«; Rq. At several H above z» ablation becomes negligible com- 
pared with sublimation so use of analytic expression ( fTSl l for n{r) 
at all heights, results in negligible error in the total energy flux. 
Over the very narrow range of 1.005 < r//?© < 1.015 around 
r - r,, Tradif) <^ '"^ increases by 1% while Tcoiiif) falls xlO^. 
Throughout the present paper, for simplicity we approximate H 
by a height-independent mean value though H changes slowly 
with depth since the temperature and hydrogen ionization vary. 
The variation in H is around a factor of 4 over the entire range 
r, -^ Rq (and far less over the much smaller range of r which 
proves to be important) while n{r) increases by a factor 10''. A 
detailed model of n(r) over a large range of r was developed by 
Porter (2007, 2008) and should be used in more precise analysis. 
In addition chromospheric n structure varies on small horizon- 
tal scales about the mean stratified form n(r) so that fluctuating 
mass loss rate is expected along the path. 



Brown et al. : Sun grazing and Impacting Comets 



From the above discussion, the total energy flux Tir) is rea- 
sonably well approximated by 



r{r) = Tradir) + Tcoiiir) 

-Lq 1 3 _( 



R<.\"^ 



r-«o)/H/£l0V 



(16) 



If we define dimensionless variables x - r/RQ,m(x) - 
M{r)IMo,p = qlR&,h = HIR^, and use Eqn. ^ for Tif), 
differential Eqn. ( fTOb becomes, with (incoming) boundary con- 
dition m{r — > oa_) - 1 



1 dm 

I'^l'^ dx 



1 



+ P- 



e 



-{x-\)lh 



x{x-pyi^ '^x^iHx-pyi^ 

where the dimensionless parameters a,[i are 
Lq 0.27 0.43 



a - 



AjtZp'^I^qRqM'J' M 



fl/3 



M 



1/3 
12 



where Mn = Mo/lO^h^Mn = MJlO'^g, 



fi 



2niimpR\v\n, 



L^ 



° 3.8x10^ 



and we define an effective mass 



(17) 



(18) 



(19) 



(20) 



In Figure 1 we show the dependence of Mol Mo on p, £.. M is 
the mass of water ice which would sublimate at the same rate 
(g/s) as the actual nucleus mass M with its specific values of 
p, £.. Note also, for dimensional clarity, that M\2 and Mn in all 
equations here are dimensionless, the dimensional scaling unit 
lO'^ g being included as a factor in the (overall dimensionless) 
numerical value of a. 

In cases where the sole processes acting on the nucleus 
were sublimation and ablation (ie hydrodynamic explosion ef- 
fects were small prior to total vaporisation) the pure sublima- 
tion/ablation solution m(x) = M(r)/Mo of Eqn. ( fTTI ) is, on the 
inbound trajectory (oo_ > x > p), 



m(x) — m_(x) 

I 



x(x — pY^^ 



1 ""i ^ ■ - 






-<-'-'^''dA 



x^'Hx-py/^ 



/2(^_^)l/2 



(21) 



while on the outbound trajectory (q < x < oo_^_) it is (with 

/rcc fX 

..dx — \ ..dx + \ ..dx) 
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Fig. 1. Conversion factor M„/M„ from ice equivalent mass M,, 
to real mass Mo versus p for X = 0.1,0.2,0.5, 1.0- Eqn. (l20l i. For 
our canonical values p - 0.5, £,- I, the true mass Mg - 4-x Mo 
since comet nuclei are less dense than ice so easier to evaporate. 



frequently in the solution. Eqns. ( I21I22| | for m(x) allow deriva- 
tion of any property of M{r). For example the mass loss rate 
M = M„m'(x).x where i = /-/Rq = (vQ/RQ)(l/p/xy^^/x^'^ by 
Eqns. (HJ. Very near the sun x —> I and is not a convenient vari- 
able since 1-p/x — > Q,m'(x) -^ oo and rates should be expressed 
via the 6 coordinate with 6 = {vq/Rq)(p^^^/x^). 

From the discussion above it is apparent that comets with 
q-Ro > O.OIRq never enter the ablation regime so can be treated 
neglecting the ablation terms. On the other hand comets with 
q - Rq < O.Ol/?0, and of large enough Mo to survive sublima- 
tion, encounter the ablation regime. If, further, they are of large 
enough Mo to survive to r = r, = LOl/?© with only a small mass 
loss fraction Am - 1 -m(xt = rt/Ro) «; 1 then their m(x) can be 
well approximated by including only the ablation terms and ne- 
glecting the preceding sublimation Am. These sublimation and 
ablation dominated regimes are treated analytically in Sections 
5 and 6. 1 while the role of ram pressure driven explosion, vital in 
the destruction of larger masses or steep entry angle, is discussed 
in Section 6.2. 



m(x) — OT+(x) 



/ p e-<-'-')/^Jx n e-^'-^^'i'dx \ 



,3 



5. Sublimation Dominated IVIass Loss Solutions 

{Mo < 10" g or q-n»H^ O.OlRa) 

5.1. Sublimation dominated solution for m(x) 

o2) 

In this case solution ( 1211 l22l i for m(x) along the incident and 

outgoing paths become 



Since Lq, Rq, Vq, Hq, H are all known, for specified p, X,, Pc, 
this solution m(x) depends solely on the the parameters q and 
Mo and we will see below that the combination Moq^^^ appears 
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Fig. 2. m{x) versus xip for various h\oP^^'^ showing cases of 
survival of sublimation through an orbit, and of terminal mass 
loss before and after perihelion. This Figure can be used for any 
values of Mo, p,p, £. because of the scaling and combination of 
variables used. 



2a 



m'f(x) = 



1 - 



2a 



as p —^ 



3p 



1/2 



(?)■" 



(23) 



which is plotted in Figure 2 in terms of m(x/p) for various val- 
ues of the parameter MoP^^^- In the following subsections we 
discuss features of this solution which are of particular interest. 

5.2. Sublimation solution properties 

5.2.1. Fractional mass lost by perihelion/after one orbit 

By Eqn. (|23]) with x - p - q/R^, the fractional masses m sur- 
viving, and Am lost by sublimation to perihelion are 



rad 



= l-A<-[l-3f^f 



0.27 



n3 



p"^M\i 



-li 



p"^M[: 



(24) 



while after one orbit they are, by Eqn. ( l23l l for m+(x — > oo) 



y,rad 
'■orbit 
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3/?'/2 
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(25) 



the last forms being for our canonical p, X. values. (These ex- 
pressions are valid for p > 7T^a^/9,4-n^a^/9 respectively). Note 
that Am - AM/M,, - AMI Mo, and that the minimum initial 
mass needed to just survive sublimation through one orbit is 8 




10 10" 

fractional mass loss A77i„ 
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Fig. 3. Value of MoP^^^ which has undergone fractional subli- 
mation mass loss Aniq by perihelion - Eqn. (|26] |. The Figure can 
be used for any values of Mo, p,p,-C because of the scaling and 
combination of variables used 

times that needed to just survive to perihehon (the correspond- 
ing minimum Og values being in the ratio 2:1). We can rewrite 
Eqn. (O as 



p'^'M, 



0.022 



[1-(1 



Am™'' 



)l/3]3 



(26) 



which lets us evaluate the initial Mo ~ 0.25Mo which will lose 
a fraction Am by perihelion - see Figure 3. 

5.2.2. Minimum initial masses and radii surviving sublimation 
to perihelion and through one orbit 

Eqns. (l24l |25]) imply the following for the minimum initial 
masses needed to just survive sublimation to perihelion and 
through one complete orbit. 
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the corresponding minimum initial sizes being 



ra^ r\ rod 

ominorhit ~ ominq 



2.7xl03cm 5.4x10^ cm 



nol/2 



■Ipp 



,1/2 



(27) 



(28) 



According to Section 2.4 for the case of a sphere this corre- 
sponds to radius bominorbu ~ aomj„orbi,/l-67 » 3.2 X lOVp'^^ 
cm. Consequently (forp = 0.5, -C = 1) any M„ < lO" g with 
p ^ I will be totally sublimated before reaching the ablation 
zone while Mg > 10" g with p < 1.01 will undergo ablation 
and be totally dominated by it or by explosion, according to the 
conditions discussed in Sections 3.4 and 6.1-6.2. 

Eqns. ( I27ll28l l can be generalised to find the mass Mo which 
is totally sublimated at any point x+ or x^ by setting m = and 
solving for Mo- Thus observation of such Xi,„d values allows nu- 
cleus mass estimation as indeed does any measurement of mass 
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loss or mass loss rate at specific x. Note that specifying m at x 
does not define M,, uniquely unless one know if it is x+ or x^. A 
small Mo has the same m at inbound point x as a corresponding 
large Mo has when outbound at x. 

From all of the above it is apparent (cf. Section 3.3) that, 
near the sun, mass loss by sublimation, and the subsequent de- 
position of kinetic energy and abundance-anomalous matter in 
the corona occur over distance scales of order Rq and timescales 
of order Rq/vq. So even total destruction of the nucleus is a slow 
extended 'fizzle', not an explosion (cf observations by Schrijver 
2011 for a case with p ~ 1.14). This contrasts sharply with low 
altitude ablation/explosion over a few scale heights H x; IQ^^Rq 
directly powered by nucleus kinetic energy - Section 6. 



5.3. Comparison with Numerical Models 

In order to see whether our simple analytic expressions above 
give reasonable approximations to mass loss, at least for q < a 
few Rq, we can compare them with those from numerical mod- 
els incorporating effects we have neglected, such as the early nu- 
merical thermal modeling by Weissman and Keiffer (1982) and 
Weissman (1983). Weissman (1983) (his Fig. 2) predicts for wa- 
ter ice spheres the change Aa resulting from sublimation through 
infall to q = R^tohe Aa ^ 2x lO"* cm for zero albedo. This can 
be compared with our expression ( l28T l with p - I namely, (with 
p = 1,X = 1 for water ice) aomin ~ 2.7 x lO-' cm. Converting 
this to the equivalent radius for the case of a sphere (Section 2.4) 
gives 1.6 X lO-' cm in good (20%) agreement with Weissman, 
given our approximations (e.g. radial insolation and no cooling 
processes). This gives us confidence in applying our analytic ap- 
proach to get estimates for the ablation regime (Section 6). 



6. Ablation and Explosion of Sun-impactors 

6.1. Ablation Dominated Solution when ram pressure can be 
neglected - Mo{l - q/r,)^'^ < M„ ^ 10'"g, r,-q>H 

For nuclei with q < r, and with original Mg large enough > 10" 
g to survive sublimation down to r < r» the residual mass 
M'g = Mo(l - Am(xt)) has (cf. Section 5.2.1) essentially the 
same velocity as if it had arrived at r» without any sublimation. 
Because of the very rapid large switch from the sublimation to 
the ablation regime, its subsequent mass loss can thus be treated 
using only the ablation term in our m{x) expression (1211 1 for a 
mass M'g falling from infinity. For simplicity here we will only 
consider Mo large enough that Am(x») «; 1 and simply equate 
M'g = Mg. Consequently our results for masses in the transition 
range around 10" g are approximate. In this Section we further 
restrict ourselves to impactors having small enough Mg and/or 
shallow enough incidence that atmospheric ablative mass loss 
is complete before ram pressure driven explosion is significant 
(exploding impactors are discussed in Section 6.2). The inbound 
evolution ofm{x) - M(r)/Mo by ablation alone can then be well 
described by using only the ablation terms in Eqn. ( |2TI ). That is 



m"'"(x) - mcoii(x) 
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The integral in Eqn. (|29] l has in general no standard closed 
analytic form but, since /j « 1 , a good analytic approximation is 
readily obtained for almost all p. Expression (|29] | is only physi- 
cally relevant down to where ablation is complete - viz J ..dx < 
3 /a/3 - which always occurs close to x = x» w 1.01 since h <& I. 
Below X - Xt the integrand increases exponentially with depth 



so is maximal at the x limit of the integral. Consequently, (un- 
like the exponential numerator) the integrand denominator can 
be well approximated by setting x,(x, - pY^^ ~ (1 - p/x^Y^^ 
(Rare cases where |x» - p\ < h need a more refined treatment, 
and inclusion of both ablation and sublimation). Solution ( |29] l 
can then be integrated to give 
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from which it follows that the nucleus will be totally ablated by 
the height 
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with corresponding atmospheric density 
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(33) 



Eqn. ( I32l ) shows that, even neglecting explosion, all nuclei 
are destroyed over a few scale heights vertically. Such short 
paths suggest that the ablation path may be approximated as rec- 
tilinear, starting around the point where n s; n, at angle (/> x (p, 
to the vertical, and of nearly constant speed Vq. From Eqns. (|4]l 
we find that 



coS(/>, = Vr(r»)/v(r,) - (1 - p/xt) 



1/2 



(34) 



If one re-derives and solves the equation for dm/dx for a 
rectilinear path one aiTives at exactly the same expressions as 
found in Eqns. ( l32l[33T l given the relation (l34l i of 0» to p/x, so 
our algebraic approximation there is equivalent to the geometric 
one of a rectilinear path. 

From Eqn. ([8]) it then follows that ram pressure driven explo- 
sion is small throughout the total ablation phase (n 
2.5 X IQ^'^Pd) cm -^ only for masses smaller than 
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10"'P^ sec^ ( 



(35) 



where we have used H% ^ 0.5 and X = l,p = 0.5. Here 
M„ is the maximum vertical entry mass which is totally ab- 
lated by depth n = n„. These values of n^"^'l,M"^/' diff'er by 
factors of a few from the estimates in Section 3.4.2 due to the 
more accurate treatment of ablation. Since x, = 1.01, for the 
range 1 < p < 1 .0099, the final sec^ </>, factor is in the range 
lO"* - 10^ and the mass above which explosion becomes dom- 
inant for such grazing impactors is in the range 3 x 10'^"'^ g 
for Pe = 1 and 3 X 10'°"'^ g for P(, = 0.1. This solar Hmit is 
much larger than for the conditions of the Shoemaker-Levy 9 
Jupiter impacts (with fragment masses estimated to be around 
10'^ g) because of Jupiter's smaller H and steep angle of those 
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Fig. 4. Estimate of the heights Ze„d and densities «(.,„/ at which 
nuclei are destroyed as a function of M^ cos"* 0, for specified f 6- 
Note that results are very sensitive to cos0, and to Pg. To the 
left of the knee ablation vaporizes the nucleus before ram pres- 
sure exceeds the strength of nucleus material. To the right ram 
pressure drives destructive explosion of the nucleus. It is clear 
that a very large range of masses is destroyed within a few scale 
heights in the chromosphere and that for all plausible masses and 
parameters, all nuclei are destroyed above the photosphere. 



impacts (~ 45 deg) which result in a much larger density n and 
hence ram pressure at a specified TV (Section 3.4.2). This mod- 
ifies somewhat the relative importance and regimes of ablation 
and explosion in nucleus destruction in the solar case, as we dis- 
cuss below after discussing the explosion regime. 



6.2. Explosive end point for nuclei of high M,, cos^^^ 0, 

For nuclei of M,, > M'^'l^f' = M,* sec^ 0,, mass loss and destruc- 
tion of the nucleus become dominated by expansion due to the 
action of the exponentially growing ram pressure. This is com- 
pounded by breakup of the nucleus and flow via hydrodynamic 
instabilities - see Chyba et al (1993), Field and Ferrara (1994), 
MacLow and Zahnle (1994), Zahnle and MacLow (1994) and 
Carlson et al (1997). Modeling this intense localised explosion 
and the subsequent fireball expansion in the solar case is be- 
yond the scope of the present exploratory paper and should be 
addressed using modified versions of the numerical simulations 
conducted for planetary cases. Here we simply make an esti- 
mate of how much the end depth is reduced by explosion and by 
the enhanced deceleration and ablation arising from the growing 
area, based on the analytic approximation given by the above au- 
thors. The explosive lateral expansion of nuclei can be approxi- 
mated as incompressible longitudinal shortening combined with 
lateral expansion with lateral size aj^ increasing as (MacLow and 
zahnle 1994 Equation (2)) 
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where a'^ is the size when pure ablation ends at A^ 
by (Equation ( l30l l) 
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with N^fiMo) (Equation (l33]l) the depth (hypothetical for M > 
M'^i") at which ablation would be complete in the absence of 
explosion. Determination of the final destruction depth in the 
explosion regime requires numerical simulation but as an upper 
limit we can evaluate the depth at which ablation across the ex- 
ploding area alone would destroy the nucleus. The relevant rate 
is 
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Combining Equations ( |33] | and ( l40b we show in Figure 4 the 
estimated density Ug,^ - Nend/H and height z„,rf = //ln(no/n) at 
which destruction ends as a function of Mo(cos 0/0.01)^ for val- 
ues of f 6 = 0. 1 , 1 . Pea is also shown to make clear how explosion 
shortens the path relative to ablation alone. The x-axis values are 
the actual Mo values for cos </> = 0.01 which maximises the abla- 
tion regime. For photospheric p - I, cos = 0.1 the mass values 
are 1000 times smaller. Thus as one decreases P(, and as p/x, be- 
comes progressively smaller than unity, the mass regime of ab- 
lation dominated destruction becomes smaller and smaller. We 
conclude that for grazing solar impacts, for Pc > 10^ dyne/cm^ 
our ablation-dominated treatment is appropriate for masses well 
above the sublimation survival limit only when the entry angle 
is very shallow (i.e. q is very near to r, = LOl/?©). For larger Mg 
or lower Pc or lower p (less shallow entry), ram pressure driven 
explosion of the nucleus becomes wholly dominant. For solid 
objects with P(, :s> I the ablation regime is much larger 

6.3. Energy deposition profiles of Sun-impacting comets 

To calculate observable diagnostic signatures of the destruction 
of comets near the sun requires modeling the primary energy de- 
position profile then calculating the thermal, hydrodynamic, and 
radiative response of the solar atmosphere - the evolution of the 
fireball which follows the primary airburst. For the impact case, 
doing this numerically along somewhat similar lines to planetary 
impact work by e.g. Carlson et al. (1995, 1997) and other work 
cited in Sections 3.4.1 and 6.2 should be done but is beyond the 
scope of the present paper. Here we simply take a first step by 
estimating the height distribution of the nucleus kinetic energy 
deposition for the ablation dominated case, touching briefly on 
the explosion domain. 

The height distribution d&jdr (erg/cm) of kinetic energy de- 
posited along the ablation path can be approximated by 



d& 

dr 



iahl 



1 dm 
* -v\r)Mo — 

2 dr 



MoV% dm Mo\% dm 
2RqX dx 2Rq dx 



(41) 
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since x is very close to unity and v very close to v© over the few 
scale heights involved. 

Using Eqn. (l30l l to find dm/dx leads to the conclusion that 
d&ldr is very sharply peaked in height z (Az ^ H) with a maxi- 
mum of 

I — I ^ 2 X 10^''(Mo/10^2)sec,;6» erg/cm (42) 

\dri,„ax 
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at 



(cm^O « 2 X 10'^(1 - plx,yi\MJW'-) 



1/2^ 



xl2s,l/3 



(43) 



In the rare cases of steep entry/higher masses we have seen 
in Section 6.2 that the onset of explosion causes rapid deceler- 
ation and disintegration of the nucleus. This will result in even 
higher volumetric heating rates than in the case of pure abla- 
tion, though numerical simulations are needed to quantify this. 
Such atmospheric power input is similar to that in solar magnetic 
flares (lO^^""'" erg/sec over loop lengths ~ 10^ cm). 

7. Discussion and Conclusions 

7.1. Summary of Nucleus Destruction Regimes 

We have shown there to be three regimes of nucleus destruction, 
depending on the values of M„,q for given values of p,X. (for 
which we use here p - 0.5, X.- I), namely 

- (a) Sublimation dominated (Section 5.2.2, Eqn. ( l27l i) 
Nuclei of p > x^,Mo < M„(p) given by 



M 10' 

M (n) - M''"" — ~ 



ominq - 3/2 



3/2 



(44) 



are completely sublimated before they reach perihelion, (or 
r»). This is by far the commonest regime for group comets. 
(b) Ablation dominated (Sections 3.4.2, 6.1) 
Since n - n„ - 2.5 x 10''*f 6 at Xn, where atmospheric ram 
pressure starts to exceeds nucleus strength, nuclei with x„ < 
p < Xt and Mi(p) < M,, < Mzip) (i.e. of sufliciently low 
mass or shallow entry angle) lose all their mass to ablation 
before explosion with 



M2{p) = MZ'(P) = 



M„Pi{i-p/x*y 



-3/2 



io'°pl(i-p/x*y 



3/2, 



(45) 



Since most larger sun-grazers have p > Xt, these are rarer 
than comets in the sublimation regime. 
- (c) Explosion dominated (Section 6.2) 

Nuclei with p < x,* and Mo > M2{p) (ie of sufficiently high 
mass or steep entry) undergo ram pressure driven explosion. 
These are also rare compared to sublimative destruction. 

These domains in the plane of (p/x,. Mo) are summarized Figure 
5 with the Mo axis scaled with respect f g. ^e = 1 represents rela- 
tively strongly bound nuclei. For fc ~ 10^ dyne/cm^ the ablation 
regime becomes confined to about Mo > 10"g, 0.99 < p/xt < 1. 
Though this is very narrow in p space, how common ablation 
domination is compared to explosion domination depends on the 
distribution of p/x,» and is very sensitive to the value of P^.. For 
very weakly bound bodies (Pg «c 0. 1 ) ablation would only dom- 
inate briefly prior to explosion totally dominating. Solid stony or 
iron bodies with P(, > 10 would, unless extremely massive, be 
wholly in the ablative destruction regime. 
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Fig. 5. Domains of destruction of comet nuclei in the plane 
(MolP\,plx^ - sin(^») for entry angle 0, = arcsin(/9/x») to the 
vertical at r = r,. NOTE the change from linear to log scale on 
the p/x, axis across the point /?/x, = 1 so as to display the very 
rapid variation in the range p/x, < 1 The axes are scaled with re- 
spect to f 6 as discussed in section 2.2. For Pe = 0.1 the ablation 
regime becomes a slender strip at 0.99 < pjx^ < 1 



7.2. Visibility of Sun-grazer destruction by sublimation 

It was shown in Section 3.3 that destruction by sublimation of 
low mass nuclei even very near the sun occurs over scales ~ Rq 
in space and ~ lO^""' s in time, making them hard to observe 
against solar atmospheric emission and scattering. However, the 
abrupt collisional stopping of sublimated material in the wake 
could raise its initial temperature by up to around 10^ K which 
may make it visible in the XUV by charge exchange or ther- 
mal line emission if it does not cool too quickly by radiation, 
conduction or expansion and if its emission measure J n^dV is 
high enough. (The first observation, in XUV, of such an event 
was made by Schrijver et al., 2011, using SDO AIA during the 
revision stage of this paper). A 10'^ g mass totally sublimated 
over a Q.SRq path into a conical wake of half angle say 10"^ ra- 
dians would have an emission measure of order 3 x lO'*'* cm"-' 
and a mean density of ~ 5 x 10^ cm"^ in a solar plasma environ- 
ment of roughly similar density. 

7.3. Signatures of Sun-impactor Destruction 

7.3.1 . Mass, Energy and Temperature of the Fireball 

The kinetic energy initially released in the lower solar atmo- 
sphere by sun impacting cometary nuclei with Mo = 10" - 10'^ 
g is 5 X (10^^ - lO"'"') erg, similar to the range of the magnetic 
energy release from 'nano' -flares to 'super' -flares. If 50% of fit,-,, 
went into heating the comet material the temperature attainable 
would be r > nipV^jAk x lO^K where k is Boltzmann's con- 
stant. This is of order the solar escape temperature since it is 
generated by infall of matter effectively from oo. 

Flare masses ejected are of the order of 10'^"'^ g, also com- 
parable with those of larger impacting cometary nuclei. Thus a 
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comet impact in the dense chromosphere should produce phe- 
nomena somewhat resembHng a solar (magnetic) flare, though 
the initial comet volume is much smaller, its density much 
higher, and the heating rise time faster {< 10 s). Because of 
the initially small size scale and high density and pressure we 
might expect rapid initial cooling (seconds) by radiation, expan- 
sion and conduction, followed by a radiative/conductive decay 
phase more like that of solar flares (~ 10^"^ s) as the fireball 
expands outward and upward, sweeping solar plasma with it (cf. 
the Carlson et al. 1997 study of the Jupiter impact fireball). 

7.3.2. Cometary Flare Size, Emission IVIeasure and Optical 
Depth 

If optically thin, the radiative output of a hot plasma of volume 
V depends on its temperature and its emission measure EM - 
Jyn^dV ~ a^iMoliicnipf ^ 10'***(M„/10'2)2/(fl/10'^)VjUc cm-^ 
for a nucleus of mass Mo when it has expanded to size a. (Here 
^c is the mean comet mass per particle in units of nip). Explosion 
of a nucleus occurs over a few scale heights » 10^ cm vertically, 
traversed in ~ 10 s at speed VqCos^. If the lateral expansion 
speed of the hot wake were similar then the lateral dimension of 
the exploding wake would also be a few scale heights. Its optical 
depth for a process with absorption cross-section cr would then 
be T a; crMJnipa^ = (MJ10^^)(o-/lQ~^°)/fic/(a/lQy. Thus 
when a nucleus of, say, M„ = 10^^ g (Oo ~ 10"* cm) becomes 
optically thin {a - 10*" cm) for cr — 10"^" cm^, it would have 
an EM « 10"*^ cm"^ for that emission. This EM is only for the 
cometary material itself, but the EM of the heated atmosphere 
may also be important since an atmospheric mass > Mg is in- 
volved in decelerating it. Since H ~\ arcsec we may expect to 
see initially a source of a few arcsec expanding at about 0. 1 arc- 
sec per sec at chromospheric altitudes and atmospheric density 
^ 10'^"^^ cm"^ depending on M,, as per Figure 4. XUV line 
spectra should exhibit highly non-solar abundances (eg 0:H ra- 
tio) as should any matter ejected into space (cf. Iseli et al 2002). 
There may also be a variety of nonthermal radio and other signa- 
tures (e.g. charge exchange - e.g. Lisse et al. 1996 - and impact 
polarized spectrum lines - e.g. Fletcher and Brown 1992, 1995) 
arising from plasma phenomena driven by the highly supersonic 
expansion, such as shock and/or turbulent acceleration of elec- 
trons and ions and charge separation currents as the ablated mat- 
ter decelerates in the atmosphere - cf. review by Pegaurier (2007) 
of fusion pellet ablation physics. 



flares and CMEs their energy density is so high that local disrup- 
tion of magnetic fields and triggering of larger scale events are 
not impossible. 

7.4. Conclusions 

Our simple analytic treatment of comet nucleus sublimation 
gives results for mass loss in reasonable agreement with previ- 
ous numerical simulations . We have proved that nuclei reaching 
p - q/Re < r,/RQ w 1.01 undergo intense bombardment by 
solar atmospheric ions (energy 2 keV per nucleon) the rate in- 
creasing exponentially with depth on a scale of around 500 km. 
It follows that, above depth r« where n - n^ - 2.5 x lO'Vcm"', 
evaporative destruction is essentially by sublimation only. On 
the other hand, the behavior of any nucleus surviving sublima- 
tion to below this height is dominated by atmospheric collisional 
eff'ects - ablation and ram pressure driven explosion. This creates 
an exploding air burst followed by a fireball spreading and ris- 
ing through the atmosphere similarly to the Shoemaker Levy 9 
impacts on Jupiter (e.g. Carlson et al. 1997). 

Because of the exponential distribution of atmospheric den- 
sity and the small atmospheric scale height H (compared to Rq) 
the terminal heights of these cometary flares is only logarithmi- 
cally sensitive to the incoming cometary mass (and to the den- 
sity and latent heat of the nucleus). It only ranges from the upper 
chromosphere for small to moderate comet masses with typically 
shallow entry angles to the low chromosphere for very massive 
steep entry comets. 

The ablated matter (atoms and small particles) decelerates 
very quickly in collisions with the atmosphere, depositing its 
kinetic energy and momentum over a few scale heights in 
tens of seconds, heating the debris and local solar plasma to 
X-ray emitting temperatures similar to solar flares and with 
comparable plasma masses and emission measures. Observing 
such very small impulsive but very luminous 'cometary flares' 
poses a fascinating challenge but will enable new diagnostics of 
cometary element abundances from the highly non-solar flash 
emission spectra. In addition the downward momentum impulse 
should generate cometary sunquakes as observable ripples in the 
photosphere akin to those found in conventional magnetically- 
energized solar flares. In the case of impacts by the most massive 
comets (10^'^ g or so) the cometary flare energy release (2 x lO-'^ 
erg) is much larger than that of the largest solar flares ever ob- 
served. An impact of this magnitude would have very significant 
terrestrial efifects. 



7.3.3. Cometary Sunquakes and other solar disturbances 

The phenomenon of flare induced sunquakes - waves in the pho- 
tosphere - discovered by Kosovichev and Zharkova (1998) and 
now widely studied (e.g. Kosovichev 2006) should also result 
from the momentum impulse delivered by a cometary impact. 
All such impacts, however small the comet mass, involve a huge 
kinetic energy density (ram pressure) pv|/2 » 10^^ erg/cm-'. 
This is ~ 10^" times the thermal energy density (pressure) of the 
photosphere or the magnetic energy density in a sunspot, being 
the energy density of a 40 MG magnetic field! Even when the 
kinetic energy is converted by ablation and ram pressure to heat 
and kinetic energy of explosion, it is initially spread over only a 
few scale heights. Hence the explosive airburst energy density of 
a comet like a Shoemaker-Levy 9 (10'^ g) is about lO'' erg/cm-' 
equivalent to that of a 5 kG field and still more than the thermal 
energy density of the photosphere. Thus, while the total energy 
of most impacting sun-grazers is small compared to that of large 
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